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1 Introduction 

1.1 The electronic many-particle problem 

The fundamental Hamiltonian for the electrons in solid-state theory has the well-known form 



In these lecture notes, we assume that the potential V{r), generated by the atomic nuclei, is 
perfectly lattice periodic. The operators ipi^\r) annihilate (create) an electron at the real-space 
position r with spin s =t, i- Despite the fact that the Hamiltonian ([T]) only describes the 
electronic degrees of freedom, the calculation of the electrons' properties poses a difficult many- 
particle problem which cannot be solved in general. The strategies to deal with the many-body 
problem can be grouped into two main categories. 

i) Model-system approaches: 

In order to explain experiments, it is often sufficient to take into account only a limited 
number of degrees of freedom in the Hamiltonian ©. Therefore, the full problem is 
replaced by a simpler 'model Hamiltonian' which describes certain electronic properties 
of a material. Celebrated examples are the Heisenberg model for magnetic insulators and 
the BCS Hamiltonian for superconductors. In many-particle theory, in general, and for 
transition metals and their compounds, in particular, multi-band Hubbard Hamiltonians 
provide the standard models, see Sec. II. 31 

ii) Ab-initio approaches: 

In order to cope with the full Hamiltonian ([T]), one has to resort to approximations which 
are necessarily cruder than those designed for the investigation of much simpler model 
systems. The most frequently used ab-initio approach is the 'Local-Density Approxima- 
tion' (LDA) to 'Density-Functional Theory' (DFT), see SecO 

1.2 Density-Functional Theory (DFT) 

The 'Density-Functional Theory' (DFT) is based on a theorem by Hohenberg and Kohn |[T1. It 
states that there exists a universal functional W{n{r)} of the electronic density n{r) such that 



has its minimum, £"0 = E{no{r)}, at the exact ground-state density no(r) of the Hamilto- 
nian (dl) and Eq is the corresponding ground-state energy. Since it is impossible to determine 




(1) 





(2) 
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the functional W{n{r)} for many-particle systems exactly, it is necessary to develop reasonable 
approximations for it. Usually, one writes 

W{n{r)} = T{n{r)} + ^ J d'r j d\' ^^j^^^ + E^,{n{r)} , (3) 

where T{n{r)} is the 'kinetic-energy functional' and _E'xc{'^(^)} is the 'exchange-correlation 
functional' which contains all Coulomb-energy contributions apart from the Hartree term that 
was separated in ([3]). Both functionals are unknown. Approximate expressions for T{n{r)} 
and E^dn^r)} are usually derived by considering a free electron gas. The kinetic energy of 
such electrons in the Hartree-Fock approximation is ~ n^^'^ where n is the density of the ho- 
mogeneous electron gas. Therefore, a common approximation for the kinetic-energy functional 
in ©is 



T{n(r)} 



10m 
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In the same way, one may approximate the exchange-correlation potential in the form 

3e2 



E.c{n{r)} = - 



dV— (Stt 



(4) 



(5) 



To work with the energy functionals & and Q is a simple example of a 'Local-Density Ap- 
proximation' (LDA) because only the local density appears in W{n{r)} and corrections, e.g., 
involving gradients Vn(r), are absent. 

The DFT in LDA, as introduced so far, provides an approximate way to determine the ground- 
state energy and the electronic density in the ground state. These quantities are of interest if 
one aims, e.g., to determine the ground-state lattice structure or lattice parameters of a material. 
Most practical applications of the DFT, however, rely on an additional concept, the 'Kohn- 
Sham scheme'. This scheme is based on the assumption that, for each system of interacting 
particles, there exists an effective single-particle Hamiltonian 



= E/ d'^v^l(^) 



2m 



+ V{r) 



(6) 



\r — r 



which has the same ground-state density no(r) as the original many-particle Hamiltonian ([T]). In 
general, one cannot prove rigorously that such a single-particle Hamiltonian exists; this poses 
the 'v-representability' problem. If a system is w-representable, however, the 'Kohn-Sham 
potential' in Q is given by 



d 



dhir) 



T{h{T)}-T' {n(r)} + ^., {n(r)} 



(7) 



h(r)=n{r) 



where T'{n(r)} is the minimum kinetic energy of free non-interacting particles with the density 
distribution n(r). Usually one assumes T'{n(r)} = T{ri(r)} such that both terms cancel each 
other in ([7]). 
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For our formulation of a self-consistent 'Gutzwiller Density-Functional Theory' (GDFT) in 
Sec. m it will be convenient to introduce a basis of local ('Wannier') orbitals (l)i^a{r) which 
are centered around the L lattice sites i and carry the spin-orbital index a. With this basis, the 
Hamiltonian ^ can be written as 

id a-,cr' 

Here, the 'electron transfer' or 'hopping' parameters 
depend on the particle density 

^(^) = EEC(^)<^.-,-'(^)(4sv)^o , (10) 

where \\Pq) is the ground state of ([8]), 

Hfl^o) = Eol^o) ■ (11) 

At least in principle, the self-consistent solution of the 'Kohn-Sham equations' (|8l)- (fTT)) is the 
central part of most DFT applications in solid-state physics. Note, however, that actual numer- 
ical implementations of the DFT usually do not work with Wannier functions but use atomic 
orbitals or plane waves as basis sets. 

Despite the rather drastic approximations which have led to the Kohn-Sham equations, a com- 
parison of theoretical and experimental results has revealed a remarkable agreement for a large 
number of materials. Therefore, the LDA has become the most important tool for the investi- 
gation of electronic properties in solid-state physics. There are, however, well-known problems 
with certain classes of materials. For example, band gaps in insulators or semiconductors are 
usually found to be significantly smaller in DFT than in experiment. Even bigger discrepancies 
arise for materials with strong local Coulomb interactions. These are, in particular, transition 
metals, lanthanides, and their respective compounds. Such systems have been investigated in 
the past mostly based on model systems, which we discuss in the following section. 

1.3 Multi-Band Hubbard models 

We distinguish 'localised' orbitals, a E £, and 'delocalised' orbitals, a E d, where the localised 
orbitals are those which require a more sophisticated treatment of the local Coulomb interac- 
tion than provided by the LDA. The natural starting point for such a treatment is a multi-band 
Hubbard model of the form 

i^H = iyo + $^i^^;c, (12) 

i 
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# 


Atomic eigenstate \r) 


Symmetry 


enersv Er 


1 


1 t,t) 


'A, 


U'-J 


2 


(1 t,i) + U,t))/v^ 




U'-J 


3 






U'-J 


4 


(1 t,i)-|i,t))/v^ 




U' + J 


5 


(1 ti,o)-io,n))/v^ 




U-Jc 


6 


(1 ti,o) + io,n))/v^ 


Ml 


U + Jc 



Table 1: Two-particle eigenstates with symmetry specifications and energies. 



This model contains a general two-particle interaction in the localised orbitals, and fixed hop- 
ping parameters, tl'J , and orbital energies, ef'°^ = t^f . Since the parameters t^j are usually 
derived from a DFT calculation, see Eq. they already contain the Coulomb interaction on 
a DFT level. For the localised orbitals this means that, through the on-site energies ' , the 
Coulomb interaction appears twice in the Hamiltonian li^c- We will address this so-called 
'double-counting problem' in Sec.|5l 

In the context of the Gutzwiller variational theory we need the eigenstates \r)^ and the eigen- 
values Ei^r of the Hamiltonian Hi-^c- They allow us to write Hi-^c as 

Hi-c = ^ Ei^rrhi,r , rhi^r = \r)ii{r\ ■ (15) 
r 

As a simple example, we consider a model with two degenerate Cg orbitals in a cubic environ- 
ment. In this case, we may set e'^''^ = and the local Hamiltonian then has the form 

Hi-c = f^e^fie^i + ^' ni^sn2,s' - ni^sn2,s (16) 

e s,s' s 

s 

where e = 1, 2 labels the Cg orbitals, s =t, i is the spin index and we use the convention 
f =1, I =|. For Cg orbitals, only two of the three parameters in (fT6l) are independent since 
the symmetry relations U' = U — 2J and J = Jc hold. In our model, we have four spin- 
orbital states a = (e, s) per atom, leading to a 2"^ = 16-dimensional atomic Hilbert space. 
All eigenstates \r)i of Hi-c with particle numbers 7^ 2 are simple Slater determinants of 
spin-orbital states \a) and their energies are 

Er = (iV = 0, 1) , 

Er = U + 2U'-J (iV = 3), (17) 
Er = 2U + AU' - 2 J (A^ = 4) . 

The two-particle eigenstates are slightly more complicated because some of them are linear 
combinations of Slater determinants. We introduce the basis 

\s,s') = ci,c{^,\()) , \n,o) = cyijo) , \o,n) = &lAJ^) (i^) 
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of two-particle states, which are used to set up the eigenstates of Hi-^c, see table[il The states of 
lowest energy are the three triplet states with spin 5 = 1, which belong to the representation A2 
of the cubic point-symmetry group. Finding a high-spin ground state is a simple consequence 
of Hund's first rule. Higher in energy are the two degenerate singlet states of symmetry E and 
the non-degenerate singlet state of symmetry Ai. 

The eigenstates of the local Hamiltonian Hi-c play an essential role in the formulation of the 
multi-band Gutzwiller theory in Sec. [2l Since in most applications only a finite (and not too 
large) number of localised orbitals is taken into account, these eigenstates can be readily calcu- 
lated by standard numerical techniques. The special case of a 3(i-shell in a cubic environment 
has been discussed analytically in great detail in the textbook by Sugano, Ref. |l2l. 

2 Gutzwiller wave functions 
The single-band Hubbard model 

To understand the main physical idea behind the Gutzwiller variational theory it is instructive 
to start with a consideration of the single-band Hubbard model 

^iB = X] ^i'i^\s^j,s + U^di , di = fii^tfii^^ . (19) 

In Hartree-Fock theory, one uses a variational wave function which is a one-particle product 
state 

l^o) = n^7lO)- (20) 

7 

in order to investigate many-particle Hamiltonians such as (fT9l ). It is well known, however, that 
such wave functions are insufficient for systems with medium to strong Coulomb interaction 
effects, see, e.g, our later discussion in Sec. HI It is a particular problem of a Hartree-Fock 
treatment that local charge fluctuations can only be suppressed in that approach by a spurious 
breaking of symmetries. Therefore, it usually overestimates the stability of phases with a broken 
symmetry. Hartree-Fock wave functions, however, can still be a reasonable starting point in 
order to set up more sophisticated variational wave functions. This leads us to the general class 
of 'Jastrow wave functions' JHIHl, which are defined as 

l^^j) = A l^^o) . (21) 

Here, |!^o) is again a one-particle product state and Pj is a correlation operator, which can be 
chosen in various ways in order to minimise the variational ground-state energy. The 'Gutz- 
willer wave function' (GWF) is a special Jastrow wave function with a particular choice of the 
correlation operator Pj. It was introduced by Gutzwiller |I5]-|3 in the form 

|^^)^P^|.Z^o) = n^'l'^o) (22) 



Gutzwiller Theory 



5.7 



and with the purpose to study ferromagnetism in a single-band Hubbard model. The (local) 
'Gutzwiller correlation operator' 



for each lattice site i contains a variational parameter g (with < g < 1), which allows one 
to supress local double occupancies that are energetically unfavourable for a finite Hubbard 
interaction U > 0. 

The Hilbert space of the local Hamiltonian for the one-band Hubbard model is four-dimensional 
where a local basis |/) is given by the states |0), |t), and \d) for empty, singly-occupied 
and doubly-occupied sites, respectively. By working with the occupation operator di in (|23l) . 
Gutzwiller singled out the doubly-occupied state \d). A more symmetric definition of the local 
Gutzwiller correlator (|23l) is given by 



where the operators rhij = are the projectors onto the four atomic eigenstates |/). The 

operator (|24|) contains four parameters \j instead of only one parameter g in Gutzwiller' s defini- 
tion (l23l) . It can be readily shown, however, that the operators (|23l) and (|24l) define the same sets 
of variational wave functions as long as the respective one-particle states \^o) are also treated 
as variational objects. Therefore, the wave functions, defined by (|24l) . contain more variational 
parameters than are actually needed. This surplus of parameters will turn out to be quite useful 
when we evaluate expectation values in the limit of infinite spatial dimensions. Moreover, for 
the multi-band generalisation of Gutzwiller wave functions in the following section, Eq. (|24l) is 
the most natural starting point. 

Multi-band Hubbard models 

It is pretty obvious [HI [3 how the Gutzwiller wave functions (|22l) can be generalised for the 
investigation of the multi-band Hubbard models (fT2l) . The starting point is again a (normalised) 
single-particle product state \\Po) to which we apply a Jastrow factor that is a product of local 
correlation operators. Hence, the multi-band Gutzwiller wave functions are given as 



P^^g'^^ = l-{l-g)d, 



(23) 




(24) 



1'G)=PG|'^0)=n^l^0) ' 



(25) 



where, as in (|24l) . we might work with a local correlation operator of the form 




(26) 



The variational parameters Xi-r^ allow us to optimise the occupation of each eigenstate \r) of the 
local Hamiltonian Hi-c. In multi-orbital systems, however, these states are usually degenerate 
and not uniquely defined. Moreover, it is not clear whether, in a solid, the (atomic) eigenstates 
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\r)- lead to the best variational ground state of the form (l25l) . Instead of (|26|) it may therefore 
be better to work with the general local correlation operator 

A= 5^ V„r/|r),,(r| , (27) 
ri,r' 

which contains a matrix \i-ri,r' of variational parameters. The analytical evaluation of expec- 
tation values, which we discuss in the following section, can be carried out without additional 
efforts for the general correlation operator (|T7l) . In numerical applications, however, we often 
have to restrict ourself to the simpler operator (l26l) since the number of parameters Aj^r^.r/ may 
become prohibitively large. Alternatively, one can try to identify the 'relevant' non-diagonal 
elements of \-ri,ri and take only these into account. Such strategies have been discussed in 
more detail in Ref. ifTOl . 

For systems without superconductivity, the Gutzwiller wave function should be an eigenstate of 
the total particle number operator 

N = Y^n,^,. (28) 
This requires that TV commutes with Pq^ which leads to 

Y,\,rA\r\-\n)\r)^^{r'\ = Q (29) 

r,r' 

where \r\ is the number of particles in the state \r)^. From equation (|29] ). we conclude that 
\i-r,r' can only be finite for states \r)^ , \r')^ with the same particle number. In a similar way, 
one can show that these states have to belong to the same representation of the point symmetry 
group. To study superconducting systems, one works with BCS-type one-particle wave func- 
tions for which the particle number is not conserved. In this case, the variational-parameter 
matrix \i-r,r' has to be finite also for states \r)^ , \r')^ with different particle numbers, see 
Refs. flaEl. 

To keep notations simple in this tutorial presentation, we will restrict ourself to the case of 
a diagonal and real variational-parameter matrix and do not consider superconducting states. 
Consequently, the local correlation operators are Hermitian, P} = Pi. Moreover, we work with 
a spin-orbital basis a for which the non-interacting local density matrix 

is diagonal, 

Ci-a^a' = Saynl^ . (31) 

This can always be achieved (i.e., for any \\Po)) by a proper transformation of the local basis a. 
To simplify the notations further, we will frequently drop lattice-site indices in purely local 
equations. 
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3 Gutzwiller energy functional in infinite dimensions 

The evaluation of expectation values for Gutzwiller wave functions remains a difficult many- 
particle problem even in the simplest case of a single-band Hubbard model. It has been achieved 
for this model in one dimension both for paramagnetic and for ferromagnetic states [[T2l - [T6ll . 
In the opposite limit of infinite spatial dimensions, expectation values can be evaluated for 
the general class of wave-functions (1251) . In this section, we summarise the main technical 
ideas behind this evaluation and discuss the resulting energy functional. An application of this 
functional to finite-dimensional systems is usually denoted as the 'Gutzwiller approximation' 
because, for the single-band model, Gutzwiller has derived the very same functional [|5]-d by 
means of combinatorial techniques ifTTl . 



3.1 Diagrammatic expansion 

In order to determine the expectation value 

of the Hamiltonian (fT2)) we need to evaluate the following quantities (i ^ j) 



{mm,,rWo) = UP,rn.,^rPr)T\P!) , (34) 
{^gWg) = (n^'). • (35) 



The r.h.s. of all three equations (I33l)-(l35l) can be evaluated by means of Wick's theorem because 
the wave function |!Z^o) is a single-particle product state. In this way, we can represent all 
contributions by diagrams with 'internal vertices' / (from operators P^), 'external vertices' % in 
Eq. (|34l) (or % and j in Eq. (l33l) ) and lines 

Plf = {claCv,a')^0 (36) 

which connect these vertices. This diagrammatic expansion, however, is still very complicated 
even in the limit of infinite spatial dimensions. As shown in more detail in Refs. [Ellll, it is very 
beneficial in this limit to introduce the (local) constraints 

1 = {Pfh.. (37) 

{cIA,.')^0 = {4,aP^^l,a')^0 ■ (38) 



These constraints do not restrict our total set of variational wave functions (l25l) because they 
merely exploit the fact that we have introduced more variational parameters A r than are actually 
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\ Ri;k — Rj;k\ , (40) 



Fig. 1: Two examples of double- 
(II) occupancy diagrams (I ^ i). Dia- 

gram (II) vanishes due to the con- 
straints 071) . Diagram (I) vanishes 
in the limit of infinite spatial dimensions. 

needed; see the discussion on the single-band model in Sec. |2l Note that moving the operator 

Pi relative to and ^, does not alter the whole set of constraints (|37l) . (|38l) . 

The constraints (|T7l) . (l38l) have an important consequence: Each diagram that results from (|33]) - 

(1351) is non-zero only when all of its internal vertices / are connected to other vertices by at least 

four lines. As a simple example, Fig. [T] shows two (first order) diagrams which result from (|34l) . 

While the constraints do not affect diagram (I), they ensure that diagram (II) vanishes. 

In order to have a meaningful (i.e., finite) kinetic energy per lattice site our lines have to vanish 

like ^ 

^^'i ^ ^/2^IK-ill ' ^^^^ 

on a hyper-cubic lattice in the limit of large spatial dimensions D. Here, we introduced the 
'New-York metric' 

D 
k=l 

where Ri-^k is the fc-th component of the lattice site vector Ri. Note that the number of neigh- 
bouring sites with distance | |i — j 1 1 is given by 

N\^-^^\\ = 2D\\'-^\\ . (41) 

The scaling behaviour (l39l) in infinite dimensions has significant consequences for our diagram- 
matic expansion. As an example, we consider diagram (I) in Fig.[il If we skip, for simplicity, 
any spin-orbital dependence of lines, this diagram leads to the contribution 

diagram (I) ~ ^ i^^^ ^ O (^1^ Q . (42) 

where the scaling 1/D results from equations (l39l) and ((4T)) . In general, one can show that 
in infinite dimensions a diagram vanishes if it contains an internal vertex that is connected to 
other vertices by three or more lines. The constraints (l37l) . (1381) ensure that this is the case for 
all diagrams which contain at least one internal vertex. Our arguments, so far, only apply to 
diagrams in which all internal vertices are connected to the external vertices i (or i and j). Of 
course, if we apply Wick's theorem to equations (I33l)-(l34l) we also obtain diagrams with internal 
vertices that are not connected to the external vertices. These diagrams, however, are exactly 
cancelled by the norm diagrams from Eq. (|35l) as shown in Ref. [|8l. In summary, we therefore 
end up with the simple results 

{clc^yha = {iM,AWc,yPj))^^. (43) 

{rhi;r)^a = ({Pirhr.rRd) , (44) 

\ / <I'o 
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for expectation values in the limit of infinite spatial dimensions. These expectation values and 
the constraints (|37|) . (1381) will be further analysed in the following section and lead us to the 
Gutzwiller energy functional in infinite dimensions. 

3.2 Energy functional for multi-band systems 
Notations 

We assume that the 2N (localised) spin-orbital states a are ordered in some arbitrary way, 
a = 1, . . . , 2N where is the number of localised orbitals per lattice site. In order to set 
up a proper basis of the local Hilbert space which belongs to Hi-c, we introduce the following 
notations for the 2^^ possible electronic configurations ('Slater determinants'). 

i) An atomic configuration / is characterised by the electron occupation of the orbitals, 

/ G {0; (1), . . . , (2iV); (1, 2), . . . , (2, 3), . . . (2iV - 1, 2iV); . . . ; (1, . . . , 2N)} , (45) 

where the elements in each set / = (cti, (T2, . . .) are ordered, i.e., it is cti < 0-2 < • • •• In 
general, we interpret the indices / as sets in the usual mathematical sense. For example, 
in the atomic configuration /\/' only those orbitals in / that are not in /' are occupied. 
The 'complement' / is defined as 

J=(l,...,2iV)\/. (46) 
where (1, . . . , 2N) is the state with the maximum number of 2N electrons. 

ii) A state with a specific configuration / is given as 

\I) = CU0)^l[^i\0)=dl...dl,,\0) , (47) 

where the operators cj. are in ascending order, i.e., it is cii < (T2 • • • < ^i/i and |/| is the 
number of particles in /. Products of annihilation operators, such as 

will be placed in descending order, i.e., with ai > a2 ■ ■ ■ > cr\i\. Note that we have 
introduced the operators C\ and Cj just as convenient abbreviations. They must not be 
misinterpreted as fermionic creation or annihilation operators. The 'sign function' 

f{a,I) = {lUa\ci\I) (49) 

counts whether an odd or even number of commutations is required to place a in its proper 
position in / (/(cr, /) = =i=l). It vanishes if cr G /. 
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iv) The operator m/ // = |/) (/'| describes the transfer between configurations /' and /. It 
can be written as 

rhij, = C\C,, n (1 - ^-") (50) 



where J = IVM' . K special case, which derives from (pOl) . is the occupation operator 

mi=\l){I\ = Wn,W{\-n,,) . (51) 

o-e/ o-'e/ 

The states |/) form a basis of the atomic Hamitonian's Hilbert space. Therefore, we can write 
the eigenstates of the local Hamiltonian (fT5l) as 

|r) = 5^T,,r|/) (52) 

with coefficients Ti^r- For a simple example, see the two-particle states in table [H 
Local energy 

The evaluation of the r.h.s. of Eq. (|44)) is straightforward if we use 

PmrP = Xprnr ■ (53) 

This equation gives us 

(mr)^G = ^r"^r , (54) 

where 

^o^ = ^rhrU=J2\Ti,r\'m'j , m^j = l[nll[{l - n^) . (55) 

Here we have used Eqs. (|26|) . (|3TI) . and (|5TI) . 

Hopping expectation values 

For the evaluation of (l43l) we start with 



r,r' ii,i[,i2A 

which follows from Eqs. (|26l) . (l50l) . (|52|) . Note that the second operator Pc^,P in (l43l) is just the 
conjugate of (|56l) with cr replaced by a'. Hence, the only remaining expectation values which 
we need to evaluate in (l43l) have the form 

E{I, J, J') = {mi.ij:mj.j^j>)q,^ . (57) 



An application of Wick's Theorem to (1571) leads, in general, to a number of diagrams with 
(potentially) several lines connecting the lattice sites i and j. At this point, however, we again 
apply the infinite-dimensional rule that all diagrams with three or more lines connecting i and j 
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can be discarded. Hence, the only remaining diagrams are those with exactly one line between 
i and j. Together with Eq. (|3TI) . we therefore find 

E{I, I'; J, J') = ^ /(7, /')5/,/'u.T^ E /(^'' W,i'T^(47S-,7')^o • (58) 

Altogether, we obtain the following result for the hopping expectation value (|43l) in infinite 
dimensions 

(cUs-,-')^G = Yl (CS'.y )^c) (59) 

7,7' 

with the 'renormalisation matrix' 

q2 = ^^$^ArAr'5^/(a,/)/(7,/0T;u.,rT^,,r'T;,^,T,,u,,r^^^ (60) 

'y r,r' i,r 

= ;^E^^V(r|ct|r)((|r)(r|c,)) . (6i) 

7 rr' " 



Constraints 

The explicit form of the correlation operator (|26l) . together with Eq. (|3TI) . gives us directly the 
explicit form of the constraints (iJTl) . (1381) 

r / 

^<x,.m° = J^A^ ^ /(a,J\a)/K,/V)T,\.,rT,V,rm?. (63) 
Summary: Structure of the energy functional 

In summary, we obtain the following Gutzwiller energy functional for the multi-band Hubbard 
models (fT2)) in infinite dimensions 

+L ^ ErXlml (64) 
r 

where, for the delocalised orbitals, the renormalisation factors are = 5^,7- The single-particle 
state enters (|64|) solely through the non-interacting density matrix p with the elements 

P(ia),(ja') = {c]^^,Ci,a)<^0 ■ (65) 

Hence, the Gutzwiller energy functional simplifies to 

(p, Ar) = E E ^f^y * PU^'U^-) + E + ^ E • 

(66) 
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It has to be minimised with respect to all elements of p and the variational parameters \r 
obeying the constraints (|62|) . (l63l) and 

= P . (67) 

The latter constraint ensures that p belongs to a single-particle product state. 
There are several ways, how the constraints (l62l) . (|63l) may be implemented in numerical calcu- 
lations flOl. In this tutorial introduction, we will simply assume that Eqs. (|62|) . (l63l) are solved 
by expressing some of the parameters Ar by the remaining 'independent' parameters A^. Equa- 
tion (|67]) is then the only remaining constraint in the minimisation of the resulting energy func- 
tion E'^^ (p, A}.) . If it is implemented by means of Lagrange parameters, see Appendix lAl the 
minimisation with respect to p leads to the effective single-particle Hamiltonian 

^o"^ = E E ^^'^1 [<) * ^lA^' + E <"'4.^.<. + E '^AA<r (68) 

which gives us \^o) as the ground state of (|68l) . 

Hfl^o) = E^\%) . (69) 
The 'fields' ri„ for the localised orbitals in (1681) are given as ifTSl 

ri. = ^E^^ (p, A^) . (70) 



The remaining numerical problem is the solution of Eqs. (l68l)-(r70l) together with the minimisa- 
tion condition 

^E--(p,A^)=0. (71) 

Numerical strategies to solve these equations have been discussed in detail in Ref. ifTOl to which 
we refer the interested reader. 

Up to this point, the effective single-particle Hamiltonian (1681) . 

Hf = Y.^^ArKr^ (72) 
k,r 

and its eigenvalues ('band-energies') E'k.r in momentum space are just auxiliary objects which 
determine |!^o) • One can readily show, however, that the non-interacting Fermi-surfaces, defined 
by the Fermi energy Ep, 

E^,r -Ef = 0, (73) 
are equal to the correlated Fermi surfaces because the momentum distribution 

^k,r = {hlAr)'^G (74) 



has step discontinuities exactly at the momenta given by Eq. (1731) . The Fermi surface de- 
fined by (l73l) may therefore be compared to those, e.g., from de-Haas-van-Alphen experi- 
ments. Moreover, within a Landau-Fermi-liquid theory, the eigenvalues E'k.r turn out as the 
quasi-particle excitation energies which can be measured, e.g., in 'angle-resolved photoemis- 
sion spectroscopy' (ARPES) experiments, see Refs. 119111911. 
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3.3 Example: single-band Hubbard model 

As a simple example we use the general results derived in Sec. 13.21 to recover the well-known 
Gutzwiller energy functional for the single-band Hubbard model Q. For this model, the (local) 
Gutzwiller correlation operator (|24l) had the form 

P = A0m0 + A|mt + A^m^ + Xad , (75) 

where 

rhdj = {1 - n^){l - hi) = 1 - - hi + d , (76) 

= hs{l - hg) = - d , (77) 

f =1, I =t, and d has been defined in Eq. (fT9l) . Equation (l54l) gives us the expectation value 
for the occupation of the four local eigenstates, 

m0 = {rhu,)^^ = - - nj) , (78) 

m, = (m,)^^ = AX(l-^°), (79) 

d = (^Vg = -^d^t^S ■ (80) 

With these equations, we can replace the original variational parameters A0, A^, A^ by their 
corresponding expectation values m0, m^, d. This simplifies the expressions for the con- 
straints (|62l) . (|63l) which then read 

1 = nifi, + m-^ + mi + d , (81) 
n° = nis + d. (82) 

Note that the second constraint (f82l) simply ensures that the correlated and the uncorrected 
(spin-dependent) particle numbers are equal, 

(^s)^G = + = = (fis)^^ . (83) 

Equations (fST]) . (|82l) can be readily solved, e.g., by expressing m0, as functions of d, 

7710 = 1 _ _ r^o + ^ (84) 

m, = n^-d. (85) 

In this way, the only remaining variational parameter is the average number of doubly-occupied 
lattice sites d. 

Finally, we can evaluate the hopping renormalisation factors (|6TI) . 

qi{d) = 5,,,.(ArfA,-n° + A,A0(l-n°)) (86) 
= '^s,^'— f====(\AM + \/w«0) , (87) 

where, in the second line, we have used Eqs. (ITSll-dSOl). In summary, we obtain the variational 
energy functional 

E«^(rf,^o) = $^(g:(rf))2 J]t„(ctc,.,)^„ +L[/d (88) 
for the single-band Hubbard model (fT9l) in the Gutzwiller approximation. 
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4 Applications 

4.1 Ferromagnetism in a Two-Band Hubbard Model 

Since Gutzwiller's ground-breaking work we know that the single-band energy functional ([88]) 
leads to ferromagnetic ground states only under very special circumstances, e.g., if the density 
of states has a sharp peak at the Fermi level and the Coulomb interaction is much larger than 
the band width. From this observation, we can already conclude that ferromagnetism, as it 
naturally appears in transition metals, is most probably related to the orbital degeneracy of the 
partially filled 3d-shell in these systems. Therefore, it is quite instructive to study ferromagnetic 
instabilities in a system with two orbitals, as a first step from the simple one-band model towards 
a realistic description of materials with partially filled 3d-shells. 

A) Model Specification 

We consider a Hubbard model with two degenerate Cg orbitals per site on a simple three- 
dimensional cubic lattice. The local (atomic) Hamiltonian for this system is given in equa- 
tion (fT6l) . We include realistic hopping parameters for transition metal energy bands to the 
nearest and second-nearest neighbours in (fT3l) . This choice avoids pathological features in the 
energy bands, such as perfect nesting at half band filling. The single-particle part of the Hamil- 
tonian (fT3l) is easily diagonalised in momentum space and leads to a density of states Do(e) that 
is shown as a function of the band filling in Fig. [21 

The case = n° = 0.29 corresponds to a maximum in the density of states at the Fermi 
energy. For this band filling, we expect the strongest tendency towards ferromagnetism. 

B) Variational Energy Functional 

For a two-band model, it is still possible to give a manageable explicit expression of the energy 
as a function of the variational parameters. The eigenstates of the two-particle spectrum all 
belong to different representations of the point symmetry group, see table HI Therefore, one 
can safely assume that the variational-parameter matrix Ar,r' = Sr,r' is diagonal and we have 
Xr = Ar' for all states \r), \r') that are degenerate due to the cubic symmetry. Then, we are 
left with 11 (independent) variational parameters uir = Xprup: 

i) two parameters for an empty and a fully occupied site: m0, /; 

ii) four parameters for singly and triply occupied sites: and tg with s =t, |; 

iii) five parameters for doubly-occupied sites: rfj'^, df^, (for the triplet ^^2), ds (for the 
doublet ^E), and dA (for the singlet Mi). 

For our degenerate two-band model, the uncorrelated local density matrix (|3TI) is automatically 
diagonal and orbital independent, 




(89) 
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Fig. 2: Model density of states at the 
Fermi energy as a function of the or- 
bital filling n„. The dashed lines indi- 
cate the half -filled case (n^ = 0.5 and 
the fillings used in this section (n„ = 
0.29 and = 0.35). The total band- 
width is W = 6.8 eV. 




As in the single-band case, the constraint equations (1621) . (1631) can be solved explicitly, e.g., by 
considering the occupations 



m0 



1 - 2n° - 2nl + dj^ + dj^ + + dA + 2dE + 4:^ + At^ + 3/ 



dr + h + 2t, + f + - {dA + 2dE + 



(90) 
(91) 



as functions of the remaining nine independent parameters. The expectation value of the two- 
band Hubbard Hamiltonian is then given by 



E. 



2b 

atom 



Y,{qlfes,o + {W - J){dV + rfi^ + rf?) (92) 

s 

+2{U' + J)dE + (f/ + J)dA + {2U + AU' - 2J){^ + + /) , 



where we introduced the orbital-independent elements 

1 



n't 



{Vtl + v^) ^ ( + 2^ + 



of the diagonal renormalisation matrix and the bare band energies 



deeDoie) . 



For comparison, we will also consider the energy 



E. 



{2b) 
dens 



5^(g:)^5.,o + (f/'-J)(rfr + <) 

s 

+2U'do + 2Udc + {2U + AU' - 2J)m + + /) 



(93) 



(94) 



(95) 



of a two-band model without the terms in the second line of the atomic Hamiltonian (|T6l) since 
this is an approximation that is often made in studies on multi-band models. In this case, there 
are seven variational parameters d\^, d\^, do, dc, t^, t^, and /, which represent the occupation of 
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the configuration states |/). The probabilities for an empty site m0 and a singly-occupied site 
rus are related to the variational parameters by 



n 







rris 

m(D - ± - - 

The renormalisation factors have the form 

1 



[c/f + t-s + 2t, + / + 4 + c?o] 



1 - 2n2 - 2n? + d}^ + + 2do + 24 + + + 3/ . 



(Is 



{y/tl + A/mJ) a/4 + a/4) 



(96) 
(97) 

(98) 



C) Ground-State Properties 



The energies (1921 ) and (1951 ) have to be minimised with respect to their respective (nine or seven) 
independent variational parameters uip and the magnetisation 



M^K-n;)/2 



(99) 



for example, by means of the algorithm introduced in Ref. [fTOll . In Fig. [3] (left), the magnetisa- 
tion M is shown as a function of U for fixed J/U = 0.2 (U' /U = 0.6). The critical interaction 
for the ferromagnetic transition, f/p'""", is about a factor two larger than its value U^^ as obtained 
from the Hartree-Fock-Stoner theory. The corresponding values f/p^"*^ always lie somewhat be- 
low the values for the Gutzwiller wave function with full atomic correlations. In general, the 
relation Miip{U) > M^snaiU) > Afatom(f^) holds, i.e., for all interaction strengths, the tendency 
to ferromagnetism is the strongest within the Hartree-Fock theory and weakest for Gutzwiller 
wave functions with atomic correlations. Furthermore, the slopes of M{U) are much steeper in 
the Hartree-Fock results than in the presence of correlations. 

The properties of the ferromagnetic phase strongly depend on the spectrum of the atomic two- 
electron configurations. To further analyse this point, we have included the case of Jq = 0, 
which changes only the excited two-electron states. A shift of the curve M{U) results towards 
smaller interaction strengths; for a given magnetisation density, a smaller interaction strength is 
required as compared to the correct symmetry case J = Jc, see Fig. [3] (left). The effect is more 
pronounced when we go to the Gutzwiller wave function with pure density correlations. These 
results indicate that itinerant ferromagnetism is strongly influenced by the atomic multiplet 
spectra. 

In Fig. [3] (left/a), we chose the particle density per band to be n° = + n^)/2 = 0.29, 
right at the maximum of the density of state curve, compare Fig. [21 In this case, there are 
finite slopes of the M(U) curves at f/p, and a 'Stoner criterion' for the onset of ferromagnetism 
applies. In Fig.[3](left/b), we chose the particle density per band as = 0.35. In this case, the 
density of states at the Fermi energy Dol-Ep.t) + -Do(-Ef,4,) first increases as a function of the 
magnetisation density. Therefore, a discontinuous transition thus occurs from the paramagnet 
to the ferromagnet. 
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Fig. 3: Left: Magnetisation density per band as a function ofU for the Hartree-Fock solution 
(HF), the Gutzwiller wave function with pure density correlations (GWdens), cind the Gutzwiller 
wave function with atomic correlations (GW atom) for (a) Us = 0.29 and (b) Us = 0.35. The 
dotted line indicates the results for GWawm with Jc = 0. The local exchange interaction is 
J = 0.2U in all curves. Right: Phase diagram as a function ofU and J for the Hartree-Fock 
solution (HF) and the two Gutzwiller wave functions (GWdens, GWatom) for (a) = 0.29 and 
(b) = 0.35; PM: paramagnet, FM: ferromagnet. 



In the case of pure density correlations, a second jump in the M{U) curve is observed that 
is absent in the other two curves. As discussed in Ref. |'20'|, this jump is related to another 
feature of the density of states. In the Hartree-Fock theory, this feature is too weak to be of 
any significance in comparison to the interaction energy. When the full atomic correlations are 
taken into account, this first-order jump at a finite magnetisation density disappears due to the 
enhanced flexibility of the variational wave function. 

Another remarkable difference between the Hartree-Fock and the Gutzwiller method lies in 
the approach to ferromagnetic saturation. In the Hartree-Fock theory, the magnetisation sat- 
urates at U values about 20% to 40% above the onset of ferromagnetism at U™ . In contrast, 
in the variational approach saturation is reached at about twice the onset value, f/sat ^ 2f/F- 
However, even when the minority spin occupancies are zero and {Sf) is constant, the majority 
spin occupancies and rfj^ vary with U since the limit of zero empty sites is reached only for 
U oo. 

In Fig. [3] (right), we display the J-U phase diagram for both fillings. It shows that the Hartree- 
Fock theory always predicts a ferromagnetic instability. In contrast, the correlated electron 
approach strongly supports the idea that a substantial on-site Hund's rule exchange is required 
for the occurrence of ferromagnetism at realistic interaction strengths. For the case = 0.29, 
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Fig. 4: Size of the local spin ((5^)^) as 
a function of the interaction strength for 
J = 0.2U and band filling = 0.35 
for the Hartree-Fock theory (HF) and the 
Gutzwiller wave functions (GWdens, GWaum)- 



the differences between the phase diagrams for the two correlated electron wave functions are 
minor. Due to the large density of states at the Fermi energy, the critical interaction strengths for 
the ferromagnetic transition are comparably small, and the densities for the double occupancies 
in both correlated wave functions do not differ much. For the larger band filling = 0.35, 
i.e., away from the peak in the density of state, the values for Up are considerably larger and, in 
the atomic correlation case, the Gutzwiller wave functions can generate local spin triplets more 
easily while keeping the global paramagnetic phase. 

The magnitude of the local spin as a function of U is shown in Fig. |4l For U 00, each 
site is either singly occupied with probability 2 — 4ra° or doubly occupied (spin S = I) with 
probability 4n° — 1. Hence, 

{{Si)^) 00 = (3/4) (2 - 4n°) + 2(4n° - 1) = 5n° - 1/2 . (100) 

For the correlated wave functions, this limit is reached from above since, for U < 00, charge 
fluctuations first increase the number of spin-one sites at the expense of spin- 1/2 sites, which 
turn into empty sites. A further decrease of U will also activate the singlet double occupancies 
and higher multiple occupancies. Thus, the local spin eventually reduces below ((S'j)^)oo. On 
the contrary, the Hartree-Fock theory does not give the proper large-[/ limit for the local spin. 
Instead, the Hartree-Fock limit is given by 

{{S,)')^ = n\3 + 2n') . (101) 

The change of ((Sj)^) at Up is only a minor effect within the correlated electron approach. In 
particular, this holds true for the case of atomic correlations, where about 90% of the local spin 
saturation value is already reached in the paramagnetic state. Again, the Hartree-Fock results 
are completely different. There, the local spin sharply increases as a function of the interaction 
strength since the absence of correlations fixes 

{{Si)y{U < t/p«P) = ((5,)^)(t/ = 0). (102) 

Finally, in Fig. [5l we display the energy differences between the paramagnetic and ferromag- 
netic ground states as a function of the interaction strength for J = O.IU. For the correlated 
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Fig. 5: Condensation energy as a function 
of U for J = 0.2U for the Hartree-Fock 
theory (HF) and the Gutzwiller wave func- 
tion (GW atom) for n = 0.29 (full lines) and 
n = 0.35 (dashed lines). 
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electron case, this quantity is of the order of the Curie temperature, which is in the range of 
100 K — 1000 K in real materials. On the other hand, the Hartree-Fock theory yields small con- 
densation energies only in the range of f/ ^ 4 eV; for larger U, the condensation energy is of 
order U. Including the correlation effects within the Gutzwiller theory, we have relatively small 
condensation energies even for interaction values as large as twice the bandwidth (U ~ 10 eV). 



4.2 Antiferromagnetic order in iron-pnictide models 

Since their recent discovery, the iron-based high-Tc superconductors, e.g., LaOFeAs, have at- 
tracted tremendous attention both in theory and experiment. From a theoretical point of view, 
these systems are of particular interest because their conduction electrons are less correlated 
than those of other high-Tc superconductors. In contrast to the cuprates, the pnictides' undoped 
parent compounds are antiferromagnetic metals at low temperatures, not insulators. However, 
the electronic mass is enhanced by a factor of two which indicates that electronic correlations 
are quite substantial in the pnictides, too. 

The theoretical description of the pnictides' normal phases already turned out to be a difficult 
problem. Standard density-functional theory (DFT) grossly overestimates the size of their mag- 
netic moment in the antiferromagnetic ground state. For example, in LaOFeAs experiment 
finds a staggered moment of m = (0.4 . . . 0.8)/iB fi2Tl - (23ll whereas DFT calculations predict 
moments of m ^ 1-8/Ub, or larger [[24ll25l . For other pnictide compounds, the comparison is 
equally unfavourable. 

The electronic structure of LaOFeAs near the Fermi energy is fairly two-dimensional and the 
bands are dominantly of iron d and (partially) of arsenic p character. A complete tight-binding 
model for LaOFeAs should therefore consist of eight bands (i.e., five iron d and three arsenic 
p bands) [|26l . see Fig. [6] (left). For many-particle approaches, however, the study of such an 
eight-band model model is obviously quite challenging. Therefore, in many theoretical works 
on iron pnictides various simpler models have been proposed to study particular aspects of these 
materials. The fact that the bands near the Fermi energy are dominantly of iron d character 
suggest the study of an effective five-band model of pure iron d-bands. Such a model has been 
proposed, e.g., in Ref. [12711 . see Fig. [6] (middle). Even simpler models may be derived if one 
only aims to reproduce the Fermi surfaces of LaOFeAs. This is achieved, e.g., by the three-band 
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model in Fig. [6] (right), which was investigated in Ref. [(28l. 

In cases where a simplified model reproduces certain properties of a material correctly, there 
will often remain doubts whether this agreement is merely coincidental or an indication that 
a model indeed captures the relevant physics of a system. A big advantage of the Gutzwiller 
theory is its numerical simplicity that allows one to study even complicated multi-band models 
with modest numerical efforts. In this way, it is possible to test the quality of simplified models 
by comparing their properties with those of more realistic Hamiltonians. In this section, we 
will compare the magnetic properties of all three models, displayed in Figs. |6l This comparison 
provides an interesting example of the dangers that lie in the study of 'oversimplified' model 
systems. Based on a Gutzwiller theory calculation it has been argued in Ref. [28 J that the three- 
band model has a relatively small magnetic moment, in agreement with experiment. However, 
as we have shown in Ref. [,29J the magnetic properties of the five-band model are very different 
from experiment and from those of the three-band model. One must therefore conclude that 
both models are insufficient in describing the magnetic properties of LaOFeAs. In fact, it turns 
out that an inclusion of the arsenic p orbitals is essential, see below. 

In many theoretical studies, the following Hamiltonian for the Hubbard interaction Hi-,, in (fT4l) 
is used, 

i/W = H^'^'^ + Hf, (103) 

^dens _ ^[/(6,6)n6,,n6/,g+ ^ ^^,y(6,6')nM^f,',.', 

b,s b(^)b' s,s' 

b{jt)b' b{jt)b';s 

Here, we introduced Us^s'{b,b') = U{b,b') — 5s,s'J{b,b'), where U{b,b') and J{b,b') are the 
local Coulomb and exchange interactions. For a system of five correlated d orbitals in cubic 
environment, however, the Hamiltonian (|103l) is incomplete [2]. The full Hamiltonian reads 
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Fig. 7: Orbital densities in Gutzwiller theory 
(full symbols ) and in DMFT ( open symbols )as a 
function ofU (with U / J = 4:) for the simplified 
local Hamiltonian He = Hc^\ see Eq. U03^ . 
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hP + H^'^ where 



t;s,s' t,s 



+ h.c. . 



(104) 



Here, t and e are indices for the three orbitals with symmetries xy, xz, and yz, and the two Cg 

y"^. The parameters in (|104l) are of the same 



X 



orbitals with symmetries u = Sz"^ —r"^ and v 
order of magnitude as the exchange interactions J(6, b') and, hence, there is no a-priori reason 
to neglect them. Of all the parameters U{b, b'), J{b, b'), A{t), T{t), S{t, t'; t", e) only ten are 
independent in cubic symmetry. In a 'spherical approximation', i.e., assuming that all orbitals 
have the same radial wave-function, all parameters are determined by, e.g., the three Racah 
parameters A, B, C. We prefer to work with the orbital averages J oc Xlfe^b' ^')' U' oc 
Ylih^u ^ °f exchange and the inter-orbital Coulomb interaction. They are related to the 
intra-orbital interaction U = U{b, b) via U' = U — 2J . Due to this symmetry relation, the three 
values of U, U', and J do not determine the Racah parameters A, B,C uniquely. Therefore, we 
make use of the atomic relation C / B = A which is approximately satisfied in solids, too. In 
this way, the three Racah parameters and, consequently, all parameters in He are functions of U 
and J . This permits a meaningful comparison of our results for all three model Hamiltonians. 
In order to test the reliability of our approach we first compare our results for the partial densities 
of the five-band model with those from paramagnetic DMFT calculations. In Fig. |7]we show 
the density of electrons in each orbital as a function of U for fixed ratio U / J = 4. The full 
symbols give the Gutzwiller result for the simplified local Hamiltonian (11031) . H = Hq + H^^^; 
open symbols give the DMFT results [|30l . Obviously, the agreement between the Gutzwiller 
theory and DMFT is very good. This comes not as a surprise because both methods are derived 
in the limit of infinite spatial dimensions. 

Figure |7] shows a common feature of multi-band model systems. The local Coulomb interaction 
induces a substantial charge flow between the bands because, for the local Coulomb interaction, 
it is energetically more favourable to distribute electrons equally among the bands. However, 
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Fig. 8: Magnetic moment as a function ofU for: left: model with eight bands [26]; middle: 
model with five bands ^T7\] (Hartree-Fock and Gutzwiller theory); right: model with three 
bands / |2^/ and local Hamiltonians H^'^^^ (dotted) and Hf (dashed) . 

the bands described by Hq are extracted from a DFT calculation whose predictions for the 
Fermi surface reproduce experimental data reasonably well. Therefore, the artificial charge 
flow as seen in Fig.|7]is clearly a consequence of the double counting of Coulomb interactions. 
Since the (paramagnetic) Fermi surface found in DFT reproduces its experimentally determined 
shape, we assume that the same holds for the paramagnetic orbital densities. For each value of 
the interaction parameters we therefore choose orbital on-site energies which lead to a para- 
magnetic ground state with the same orbital densities as in DFT. Note that a more sophisticate 
calculation of orbital densities requires the self-consistent Gutzwiller DFT scheme which we 
shall introduce in the following section. 

In Figs. [8] we display the magnetic moment as a function of U for all three model systems. As 
mentioned before, the three-band model shows relatively small magnetic moments over a large 
range of Coulomb- and exchange interaction parameters, see Fig. [8] (right). This is in stark 
contrast to the results for the five-band model in Fig. [8] (middle). Here we find a transition to 
an antiferromagnet with a large moment, m > 2/iB which is as abrupt as in a corresponding 
Hartree-Fock calculation. However, if one takes into account the arsenic p bands, the magnetic 
moment is significantly smaller and in the range of experimental without the need of 'fine- 
tuning' the Coulomb- or exchange interaction parameters, see Fig. [8] (left). In summary, we can 
conclude that a proper understanding of the magnetic order in LaOFeAs requires the study of 
an eight-band model of iron d and arsenic p bands. It it possible, of course, that other properties 
of this compound, e.g., the superconducting order, may be correctly described by simpler model 
Hamiltonians. 



5 The Gutzwiller density-functional theory 
5.1 The Gutzwiller DFT equations 

The model-based Gutzwiller method, which we have used in the previous section, ignores the 
fact that the hopping parameters ^ are actually functions of the density n(r). Taking this 
functional dependence into account defines the 'Gutzwiller density functional theory' (GDFT). 
The explicit inclusion of the local Coulomb interaction within the Gutzwiller theory leads to 
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changes of the particle density n{r) for three reasons: 

i) The particle density n(r) in the Gutzwiller-correlated ground state 



iy^j o',o'',7.7' 



differs from the corresponding DFT expression (flOl) . 



ii) Due to the renormalisation factors in (|66l) there will be an energy gain from changes of 
the hopping parameters t'^'J = t^j {n{r)}, which requires the re-adjustment of n{r). 

iii) The Coulomb interaction can lead to drastic changes of the occupation numbers = 
P(ia),(ia) the localiscd orbitals, e.g., when the ground state is magnetically ordered. This 
also changes the non-local elements of the single-particle density matrix p and the particle 
density (I105D . 

These correlation-induced changes of the particle density are taken into account in the self- 
consistent GDFT by including the dependence of the hopping parameters on n{r). Equa- 
tion (11051) shows that n{r) and, consequently, tl'"^ are unique functions of p and \r- Therefore, 
the GDFT energy functional has the form 

(t,(t',7,7' i^j i,(TGd r 

(106) 



We assume again that the constraints (|62l) . (1631) are solved by expressing some of the pa- 
rameters Ar by the remaining 'independent' parameters Xp. The resulting energy functional 
^GDFT ^ minimised with respect to the density matrix p and the independent 
variational parameters A}., 

-^E^^^^{p,Xr)=^ , (p, A^) = , (107) 



with the usual constraint (|67l) for the non-interacting density matrix p. The minimisation with 
respect to p leads to 'renormalised' Kohn-Sham equations of the form ([8]), (fTTI) . and (11051) 
with Eq. dH]) replaced by 

iy^j cr,(7',7,7' i,aGi 

and by [[H 

respectively. The set of Eqs. (11051) . and (I106I) - (I109I) . which have to be solved self-consistently, 
constitute the GDFT. This approach was first proposed in Refs. [|3TI - [33l and has been applied 
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to various systems in Refs. [|32l - [36l . In all these works, the authors report a remarkably better 
agreement with experiment than it could be obtained by a model-based Gutzwiller calculation 
for the same materials. 

One of the main advantages of the DFT is its 'ab-initio' character, i.e., the absence of any 
adjustable parameters. Unfortunately, this benefit of the DFT cannot be fully maintained in 
the GDFT because that would require the calculation of the two-particle interaction parameters 
^cri,cr2,cr3,CT4 ^^YQ locallscd orbltals from first principles. The straightforward ab-initio solu- 
tion of this problem, namely to calculate these parameters from the Wannier orbitals (pi^air), 
is known to yield values which are much too large. Apparently, screening effects decrease the 
Coulomb-interaction parameters significantly. These effects, however, are not well understood 
and a quantitatively reliable technique for the calculation of screened Coulomb parameters does 
not exist. For this reason one usually applies the same strategy as in model-based calculations 
where the matrix elements f/f^i '^2,^3,^4 somehow parameterised, e.g., in spherical approxi- 
mation by means of a few Racah parameters. These are chosen to obtain the best agreement 
with experiment. In this context, it is a big advantage that the GDFT provides one with more 
data, e.g., with structural properties, that can be compared to experiment, see Sec. 15.21 
As mentioned before, the local Coulomb interaction appears twice in the Hamiltonian ([14)) be- 
cause it also affects the on-site energies ef '^ . There are several ways to overcome this 'double- 
counting problem' which have been proposed in the literature, see, e.g., Ref. Il37l . According 
to Refs. [|32l - [36l . the subtraction of the mean-field operator 



from Hi-c leads to good results within the GDFT. 

5.2 Application 

As an example for the relevance of the GDFT we show results on the iron-pnictides LaOFeAs 
and BaFe2As2 which have been presented in Ref. [|35l . The failure to describe the magnetic 
order is not the only problem the DFT faces in its calculations on iron-pnictides. There are 
also substantial deviations between the experimental results and DFT predictions on lattice 
parameters, in particular for the distance between Fe and As. Taking correlations into account 
more properly, as it is done within the GDFT, changes these lattice parameters significantly. 
Figure [9] shows the interlayer distance (i|eAs ^ function of J for several values of U. In both 
systems, the exchange interaction clearly plays an important role and needs to be included in 
order to reproduce the experimental value for (ipg^s. As a consequence, other properties are also 
changed significantly as a function of J, see Ref. Il35ll . It should be noted that the calculations 
in Ref. [|35l were carried out with the simplified local Hamiltonian H^'^^^ in Eq. (11031) . Taking 
the full atomic interaction into account may therefore change results, at least quantitatively. 
Nevertheless, these results already illustrate how important it is, in studies on transition metal 
compounds, to treat the local Coulomb interaction in a more sophisticated way then provided 
by state-of-art DFT methods. 
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Fig. 9: Interlayer distance d^^xs LaOFeAs and BaFe2As2 as a function of J for several 
values ofU. 

6 Summary and Outlook 

In this tutorial presentation, we have provided a comprehensive introduction into the Gutzwiller 
variational approach and its merger with the density functional theory. Numerically, the Gutz- 
willer method is rather 'cheap' as compared to other many-particle approaches. It will therefore, 
quite likely, emerge as an important tool for the improvement of existing ab-initio methods. 
There are two more recent developments which we shall briefly mention as an outlook: 

• The time-dependent Gutzwiller theory 

The Gutzwiller theory, as introduced in this presentation, can be used for the calculation 
of ground-state properties and of quasi-particle energies in the Fermi-liquid regime. For 
the description of experiments one often needs to calculate two-particle response func- 
tions such as the magnetic susceptibility or the optical conductivity. This is achieved 
by the so-called 'time dependent Gutzwiller theory'. This method is derived in a very 
similar way as the 'random-phase approximation' can be introduced as a time-dependent 
generalisation of the Hartree-Fock theory. It was first developed for single-band Hubbard 
models by Seibold et al. [13811395 and has been applied with astonishing success to quite a 
number of such models and response functions [|40l - [50l . Recently, the method has been 
generalised for the study of multi-band models II51[|52II . 



Beyond the Gutzwiller approximation 

As we have demonstrated in this presentation, the energy-functional which we derived in 
infinite dimensions (i.e., the Gutzwiller approximation), constitutes already a major im- 
provement over, e.g., the Hartree-Fock theory. It is well-known, however, that the limit 
of infinite spatial dimensions has some severe limitations. For example, if we consider 
the Fermi surface of a single-band Hubbard model, it will be independent of U as long as 
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no symmetry-broken phases are considered. This cannot be correct in finite dimensions 
which becomes evident already from straightforward perturbation theory in U [|53l . It is 
also known from a numerical evaluation of Gutzwiller wave functions in two dimensions 
that, for sufficiently large values of U, the variational ground states can be supercon- 
ducting [|54l[55|| . This is also not reproduced within the Gutzwiller approximation. In a 
recent work, we have therefore proposed an efficient diagrammatic method for the eval- 
uation of Gutzwiller wave functions in finite dimensions [|56l . It has enabled us to study 
correlation-induced Fermi-surface deformations [|56ll as well as superconductivity (un- 
published). The numerical efforts of this method are still moderate and the investigation 
of more complicated multi-band models will therefore be feasible in the near future. 
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Appendix 

A Minimisation of functions with respect to non-interacting 
density matrices 

We consider a general function E{p) of a non-interacting density matrix p with the elements 

P7,7' = (CyS)'2*o • (111) 

The fact that p is derived from a single-particle product wave function \<Po) is equivalent to the 
matrix equation fp = p. Hence, the minimum of E{p) in the 'space' of all non-interacting 
density matrices is determined by the condition 

T7^^(P)=0, (112) 

CP7',7 

where we introduced the 'Lagrange functional' 

L{p) ^ E{p)-Y,^iAf-p\m,i (113) 

= E{p) - ^ [2i.,n ^ Pm,pPp,l - Pm,l^ (1 14) 

l,m p 

and the matrix i? of Lagrange parameters i?i,m- The minimisation of (II 131) leads to the matrix 
equation 

H = pf2 + f2p- n (115) 
for the 'Hamilton matrix' H with the elements 

H^,Y = j^E{p). (116) 
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For density matrices which satisfy = p, Eq. (|115l) leads to [H, p] = 0. Hence, H and p must 
have the same basis of (single-particle) eigenvectors and, consequently, |^o) is the ground state 
of 

Hf = Y,H.y,yc%, . (117) 

7:7' 
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